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Abstract 

The Shubnikov - de Haas effect in quasi-two-dimensional compounds is studied. The 
conductivity is calculated from the Kubo formula. Two effects - the field-dependent phase shift 
of the beats and the slow oscillations of the conductivity are explained and calculated. The 
results are applicable to the strongly anisotropic organic metals and other layered compounds. 

1 Introduction 

The magnetic quantum oscillations (the de Haas - van Alphen (dHvA) and the Shubnikov - de Haas 
(SdH) effects) have been discovered long ago and were used a lot as a powerful tool of studying 
the geometry of Fermi surfaces and other electronic properties of different metals [jlj . Last years 
the quasi-two-dimensional (Q2D) organic metals attract a great interest and many works have 
been devoted to the study of magnetic quantum oscillations in these compounds (for a review 
see e.g. ||). The quantum oscillations of the magnetization is a thermodynamic effect that is 
completely determined by the density of states distribution. An exact calculation of the density 
of states is a very complicated problem but the scmi-phcnomenological theoretical description of 
the magnetization oscillations in Q2D compounds was recently provided in a number of papers 
@j H [| @i H H> 0- T nc cnem i ca l potential oscillations and the arbitrary electron reservoir due 
to the open sheets of the Fermi surface make no principal difficulties fl(i||. Of course, if the e-e 
m ' interaction drastically changes the ground state (like the FQHE effect) the problem becomes much 
^ , more complicated. But when the number of the occupied LLs is very large (tip > 100 as in most 
5_i ' of Q2D organic metals) the effect of the e-e interaction is reduced (as in the Fermi liquid) and 
5^ i can be taken into account via the renormalization of the effective mass. Another open question 
in the theory of the Q2D dHvA effect is the exact shape of the Landau levels. It depends on the 
particular type of a compound (on the type of impurities, their distribution, the interlayer transfer 
integral etc.) and makes only some quantitative differences. 
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The theoretical description of the Q2D quantum conductivity oscillations was not as successful 
although some works on this subject have appeared last years |ll} [l^, [l3[ There are still 
some open qualitative questions. One of them is the origin of the phase difference in the beats of 
the SdH and dHvA effects, ft takes place even when the harmonic damping is very strong and 
only the first harmonic is seen. An explanation and an approximate theoretical description of this 
phenomenon as well as the more extensive experimental study of this effect has been proposed 
recently JI3|. 

Another very interesting phenomenon is the slow oscillations of the magnetoresistance that have 
been observed in a number of Q2D organic metals §, || ||, [I?], |9[ |§ . The behavior of the 
slow oscillations resembles that of the SdH effect that has lead to a suggestion of additional, very 
small Fermi surface pockets in these materials. However, band structure calculations (basically 
giving a good description of the electron band structure and Fermi surface topology in the organic 
metals) have shown no evidence of such small pockets in any of these compounds. Moreover, 
while the slow oscillations are often very pronounced in the magnetoresistance, sometimes even 
dominating the oscillation spectrum, no analogous observation in oscillating magnetization (dHvA 
effect) has been reported thus far. Certainly, the slow oscillations must contain a useful information 
about the compounds but to extract this information one needs some theoretical explanation and, 
desirably, a quantitative description of the phenomenon. 

In this paper we calculate the conductivity starting from the Kubo formula and handling 
carefully with all oscillating quantities in the expression of the conductivity. The proposed theory 
explains both the slow oscillations and the phase shift of the beats and gives some quantitative 
description of these phenomena. It should be noted, however, that our analysis has one important 
limitation - it is developed only when the magnetic field is perpendicular to the conducting layers. 
When the magnetic field has substantial tilting angle our results give only a qualitative description. 

In section II we calculate the conductivity starting from the Kubo formula. The calculations 
are described in detail so that they can be easily followed. In Sections III and IV wc study the 
limit of strong harmonic damping in more details. In Sec. Ill we give the description of the phase 
shift of beats. In Sec. IV the slow oscillations are considered. 



2 The calculation of the conductivity 

First, we shall consider the SdH effect in the external magnetic field perpendicular to the conducting 
planes: H\\z. Then the electron spectrum of Q2D electron gas in magnetic field is given by 

e ntkz = Huj c (n+ -) - 2tcos(k z d) (I) 

where uj c = eB z /m*c is the cyclotron frequency, k z is the wave number in z-direction perpendicular 
to the conducting layers, d is the interlayer distance and t is the interlayer transfer integral. This 
energy spectrum can be easily obtained in the tight binding approximation after we choose the 
electromagnetic vector potential in the form A — (0, B z x,0). The interlayer bandwidth W is 
related to the interlayer transfer integral as W = At. We consider the case ep ^ TujJ c . Then the 
electron-electron interaction does not completely change the energy spectrum (as in the extremely 
quantum limit of FQHE regime) and we can take the energy spectrum (Q) as a zeroth approximation 
and then treat the impurity scattering using the perturbation theory. If magnetic field is tilted 
with respect to the conducting planes the formula becomes no more valid. This is the most 
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important limitation of the analysis that follows. Hence, at the substantial tilting angles our theory 
gives only qualitative predictions. 

Before diving into the calculations let us describe on the physical level the origin of the slow 
oscillations. They come from the interference of the oscillations of the different quantities, on 
which the conductivity depends. These quantities are the electron relaxation time, the mean 
square electron velocity (summed over all states on the FS) and the Dingle temperature. The 
product of two oscillating factors gives a constant term that is just one half of the product of the 
amplitudes of these oscillations: (1 + a cos x)(l + b cos x) = 1 + (a + b) cos x + (ab/2) cos 2x + ab/2 
(where in our case x = f^p)- The last term ab/2 is responsible for the slow oscillations. In our case 
the amplitudes a, b of the fast oscillations (for example, of relaxation time and the mean square 
velocity) are themselves the slowly oscillating functions of the magnetic field due to the beating. 
Hence, the conductivity contains the terms ~ R 2 D cos (47ri//iw c ) 2 = (i?|)/2)(l+cos (8irt/huj c )). The 
term {R 2 D /2) cos (%Trt/ftw c ) describes the slow oscillations with the frequency equal to the double 
beating frequency. The amplitude of the slow oscillations contains the square of the Dingle factor. 
Nevertheless, it can be larger than the amplitude of the fast SdH oscillations. The slow oscillations 
depend only on the transfer integral t and are not sensitive to the exact position of the chemical 
potential. Hence, the amplitude of the slow oscillations does not contain the temperature smearing 
factor and is larger than the amplitude of the Shubnikov oscillations at high enough temperature 
(compared to the Dingle temperature). 

The magnetization does not have such a slow oscillating term because in the lowest order 
on damping factors [see Eq. ( |22| ) below; the term (1 + 6 cos a;) comes, for example, from the 
Dingle factor] it has only the interference of a sin a; (1 + 6cos.x) = 1 + a sin a; + (ab/2) sin 2a; that 
does not have a constant term ab/2. Hence in the same lowest order on the damping factors the 
magnetization does not reveal the slow oscillations. 

The origin of the phase shift was explained recently in Jl4j basing on the Boltzmann transport 
equation. In this paper we shall give more strict calculation of this effect. 

To calculate the conductivity we shall use the Kubo formula [pl| . The procedure is similar to 
that in three-dimensional metals without magnetic field (|2~l]], § 7.1.2). In magnetic field only the 
new set of quantum numbers to = {n, k z , k y } instead of the momentum p and the different disper- 
sion relation (Q) should be used. The evaluation of the Kubo formula without vertex corrections 
gives 

^ = ^,» / ^A 2 (m,e)(-n' F (e)) (2) 

m 

where the volume V is to normalize the sum over quantum numbers to , e is the electron charge, 
the limits of the integral over e are (— oo;oo), n' F (e) is the derivative of the Fermi distribution 
function: 

-n' F (e) = l/{4Tcosh 2 [(e - n)/2T]} (3) 

and A(m, e) is the spectral function that is related to the electron Green's function G R (m, e) or to 
the retarded self energy part T, R (m, e) : 

„», , — 2Im£ fi (m, e) , . 

A(m, e) = -2lmG R to, e = o (4) 

[e - e(m) - ReT, R (m, e)f + [ImZ R (m, e)f 
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The formula (||) is close to the corresponding formula without magnetic field [|pH, formula 7.1.10] 
until the self energy part £ (m, e) is specified. It arises mainly from the impurity scattering. The 
constant part of the real part KeT, R (m, e) of the self energy does not influence the physical effects. 
It produces only a constant shift of the chemical potential. The small oscillations of Re"E R (m,e) 
(that are <C huj c can also be neglected because the final answer is not sensitive to ReE ii (m, e). 
Hence, for simplicity, in the later formulas we shall put ReS fl (rn, e) = 0. 

The imaginary part of the self energy ImS R (m, e), on the contrary, is very important since it 
describes the momentum relaxation of electrons. The main contribution to the resistivity comes 
from the short range impurity scattering. The matrix element of the point-like impurity scattering 
is the same for all transitions without change of electron energy. Hence, the scattering rate l/r(e) 
is proportional to the number of states to which an electron may scatter, that is to the density 
of states at the given energy. So, if one considers the point-like impurity scattering in the Born 
approximation, the imaginary part of the self energy is proportional to the density of states: 

-ImS R {m,e) = C x p(e) (5) 

The constant Cki CiU 2 l^ z where Ci is the impurity concentration, U = J V(r)d 3 r characterizes 
the strength of the impurity potential and I^, z is a number that depends on the electron wave 
function and the impurity distribution along the z-axis. The physical reason for the proportionality 
is clear. l/r(e) ~ /o(e). The formula (|5|) violates if one considers the next terms on the small 
parameter U Nll/HuJc = f /d, where Nll is the LL degeneracy per unit area, / - is the scattering 
amplitude (which is constant at small wave vector q <C l/^oi r o is the range of the impurity 
potential) and d is the interlayer distance. Usually this parameter / jd is really small. 

The constant C in (^J) contains many unknown parameters. We shall not calculate them here 
but note that C is simply related to the average Dingle temperature Tq: 

(\lmZ R {m,e)\) =C- (p(e)) = C ■ {N LL /fm c ) {I + n R ) = nk B T D (6) 

where (..) means the average value of something, ks — 1.38-10~ 16 erg/ K is the Boltzmann constant 
and nji is the density of reservoir states, that exist in many organic metals due to the open sheets 
of the FS. 

In the extremely 2D case (hw c ^S> t) the substantial deviations from the formula (|5|) are possible 
because the strong degeneracy of the LLs makes the simple perturbation theory not applicable. 
Since we consider the case 2i > ftw c (when the beating of the oscillations is visible) we shall take 
(^|) for our subsequent calculations that are now straightforward. 

Performing the summation over k y in (|^) and changing the integration over k z by the integration 
over energy e(n,k z ) we get 

w d(k z d) 2 f de 



*z. = ^"E 2 ^ ^vl(k z )J ^(e(k z ,n),e)(-n' F (e)) = 

= e*N LL d J 2^X>( £ »| J ^-A 2 (e\e)(~n' F (e)) (7) 

where the velocity v z (e, n) is determined from the dispersion relation (|l]) 

de(n,k y ,k z ) ltd . 

v z (e,n) = —f = —sm(k z d) = 

hok z h 



^\/4t 2 - (e- huj c (n +1/2)) 2 (8) 
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To go further we have to transform the sum over LLs to the sum over harmonics. This can be 
done using the Poisson summation formula (Appendix A). Substituting ( p2[ ) to (^) we get: 

,, T f de' , . i, td 2 (2-Kike'\ ( 4irkt\ f de ,n , , w , , 

'» = bNll 2 ^ E(-d w- p (^7j Jl (^:jy^ (e,e)( " <(e)) = 

fc= — oc 

k— — oo 

where for the zeroth harmonic k = one should use the expansion J\{x) — x/2 at i < 1, and the 
integral iz(e, /c) over e' can be easily evaluated with the spectral function (|4|): 

de'n,, . I 2 / nike l 



I z (e,k) = J ^A 2 (e', e )exp 



de' ( -2Im£ R (e) Y f2nike' 
' " exp 1 



/2wike\ ( 1 27rfc\ r> , 



2?r \[e-ef + [ImE R (e)Y I \ 
1 

y fku c ) \\tiaYj R {t)\ hhj c 
where 

R D (k,e) = exp (-2tt |lm£ R (e)| /hjj c ) (11) 

has the form, similar to that of the usual Dingle factor i?o(fc) = exp (— 2ir 2 k kBTu/hco c ) ■ Collecting 
the formulas and |l0|) we get: 

(-l) k 2td 2 ( Ankt\ 



2,r f de , , . (-l) fc 2i^ (4fKkt\ 

A;— — OO 



f2mke\ ( 1 27rfc\ „ , 

xex H^j(^M + ^:j^' e) (12) 

We have to consider all the terms in the expression (|l^) for the conductivity that make an 
essential contribution to the oscillations. Besides the directly oscillating term (exp(2irike/?iuj c ) 
these are the imaginary part of the self energy |lmE ii (e)| and the Dingle factor Rn(k,e). The 
imaginary part of the self energy is given by (||) where the density of states is 

= r E A ^ E ) = T^k E / 2 ,fn ^ «m),E) (13) 
2lT m (2tt) fr?r' J ««(e(m),») 



The sum over LLs in ( |13[ ) can be again represented as a harmonic series using the Poisson sum- 
mation formula. Substituting (Q) into ( p"3| ) we get 

^ , ,, k N LL ( Ankt\ f de(m) ( 2mke(m)\ 

^ll ^ , a. /47rfct\ [2mkE\ , 
= it E^o^jexp^j^^) (14) 
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From the formulas (||, || and |l4|), using Jo(0) = 1 we get 



|lmS R (m,e)| = nk B T D h + 2 £(-!)* J cos (j^j Ro{k,e) 



(15) 



fc=i 

Together with ( |lT|) this is a nonlinear equation on ImS R (m, e). We shall obtain explicit results in 
the next two subsections only in the limit of strong harmonic damping where an expansion over 
the harmonic damping factors in Eqs. ([ll]) and ( |l5| ) is possible. 

Substituting (|l5|) into (|l2|) we obtain the following expression for the conductivity: 

o zz = e 2 N LL / de(-n' F (e)) c -±- )—^{ 

^ 7rfc B T D (l + 2Er =1 (-l) t ./o (tjf ) cos (^) R D (k, 

9 /" , , , ,u 4£cP ,i. / 4irkt\ f2Trke\ ., . , , 

+e ^ 7 * ^ g^ 1 ) Jl (^7 J cos J e) (16) 

Generally speaking, one should also take into account the oscillations of the chemical potential, that 
are given by (JTo|, formula 5). The formula jl6| ) is quite huge and without further simplification is 
applicable only for the numerical calculation of the conductivity. In the next two sections we shall 
analyze the limit of strong harmonic damping. 

3 The first order on damping factors; the phase shift of the 
beats. 

For simplicity we now assume the harmonic damping to be strong and neglect all the terms that 
have additional powers of the damping factors: the Dingle factor R B or the temperature smearing 
factor 

R T (k) = (2kir 2 k B T/huj c ) /sinh (2kir 2 k B T/huj c ) 

that will appear after integration over energy with the Fermi distribution function. So, we rest 
only the first harmonic in the expression (|l^) for the conductivity. This approximation is sufficient 
to obtain the phase shift of the beating. In this approximation the Dingle factor R B (k,e) and the 
chemical potential can be considered as a constant. The conductivity then becomes 



e 2 N LL f , 2td 2 1 2t 1 

t zz = — de(-n F (e))— x j— -r j—t (17) 

2, J ft Tk B T D to c l-2coB(^)j (jg)jfc 



e 2 N LL f de ,_ n , (c)) 2^ 2 1 2co f^e\ Ji ( '47rt\ L { 2v 2 k B T D \ R ^ 
2ir J F ft irk B T D \hw c J \hw c ) \ frw c j 

Since we consider the limit Ro <C 1 the constant term is much larger than the oscillating term in 
the denominator in the first line of Eq. (|l7|). Hence, one can move up the small oscillating term 
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using 1/(1 + a) ks 1 — a, a <C 1. Then the integration over e can be easily performed. One gets 

RtRd 



e 2 N LL (AtdY f + 2cog /2tt^\ 



flJjJr I iln \ fVjJr I \ TilOr 



(18) 



where \i is the chemical potential and the temperature smearing factor Rt is given by the usual 
L-K expression: 

2TT 2 k B T/hw c 

Rt = ~~r 



sinh (27r 2 k B T / hu c ) 



If the transfer integral is still large enough Ant ^> hw c one can use the expansions of the Bessel 
function at large value of argument: 



J (x) ks y/2/nx cos (x - 7r/4) , a; > 1 (19) 
Ji(x) f» \/2/n;sin(i-i/4) ,i>1; 



Then the expression in square brackets of (J18J) simplifies and making use of the standard trigono- 
metric formulas one can write the oscillating part of the conductivity as 



e 2 N LL (AtdY f27r/A / huj c r—^ ( 

— — — x 2 cos - — R t Rd\ ttt- V 1 + a cos 



^ + 0^ (20) 
where 

Huj c ( 2ir 2 k B T D \ , , 

a = 1 H and = arctan (a) (21) 

2-7rf \ na; c / 

One important difference of this formula from the usual Lifshitz-Kosevich expression is the phase 
shift <f> of the beats of the oscillations. This phase shift can be experimentally observed by com- 
parison with the dHvA oscillations, that possess the phase offset consistent with the standard 
Lifshitz-Kosevich formula and are proportional to 



M ~ sin 



2nfi 



hljJr 



R D (T D )R T (T)cos(Airt/hu> c ~ir/A) (22) 



The above result concerning the phase of the beats of the SdH oscillations has been recently 
qualitatively confirmed by a comparative experimental study of the SdH and dHvA oscillations in 
an organic metal /3-(BEDT-TTF) 2 IBr 2 @- 

The main limitation of the proposed analysis is that the magnetic field is taken to be perpen- 
dicular to the conducting layers. A finite tilting angle of the magnetic field with respect to the 
normal to the conducting planes may be approximately taken into account by rescaling the Landau 
level separation: uj c — > lu c cos6> and the warping of the Fermi surfaceQ : t(6) = t(0) Jo(kFdtan9), 
where ftp is the in-plane Fermi momentum. But this is only a semiclassical approximation based on 
the assumption that the FS remains the same. Actually, the tilting of the magnetic field changes 
the dispersion relation and a more profound study of the effect of the tilting magnetic field on 
the transport properties is required. The quantum mechanical calculation of the dispersion rela- 



tion in tilted magnetic field in the first order on the transfer integral) 25 gives the result close to 
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that of Yamajip4||; the correction is that the Bessel's function should be replaced by the Laguerre 
polynomial J25| : 

J (k F dtaxid) -> e~ (A/2)2 x L n (A 2 /2) , where A 2 /2 = d 2 mw c tan 2 (0) /2h and n = n/hiu c 

But the use of the perturbation theory assumes that the transfer integral t is much less than the 
Landau level separation huj c . Here we consider the opposite case 2t > huj c where the beating can 
be seen and some more profound study of the effect of the tilted magnetic field is needed. For 
example, the effect of the tilting angle on the conductivity should not be described just by this 
replacement near the zeros of Jo^kpdt&iiO). So, our analysis is applicable to a nonzer o tilting 
angle only qualitatively. 

Other errors may come from the approximate expression for the self energy (^|). The Born 
approximation in Q2D case works quite well, but maybe other scattering mechanisms (especially 
for the calculation of the DoS) should be taken into account. An accurate study of this problem 
may depend on a particular type of a compound in hand. 

The above analysis does not take into account the vertex corrections. In our case this is valid 
because, according to the Ward identity, the vertex T(m, E) = p + m VpS fl (m, E). Hence, if the 
retarded self energy depends only on energy, the vertex corrections are zero. The fact that E fl ( m, e 
is approximately a function of only the energy e is a consequence of the short-range (or point- like) 
impurity potential. More precisely, if one takes a point-like impurity potential and neglects all 
diagrams with the intersections of the impurity lines in the self energy, then after averaging over 
randomly and uniformly distributed impurity positions one obtains that T, R (m, e) = S fl (e). The 
neglected graphs with the intersections of the impurity lines describe the coherent scattering on two 
impurities simultaneously. The contribution of such a scattering is small at large enough interlayer 
transfer integral. In the three-dimensional case without magnetic field the vertex corrections 
produce the factor (1 — cos a) instead of 1 in integrand for the transport scattering relaxation time 
where a is the scattering angle. But if the scattering probability is independent of the scattering 
angle (as for point-like impurities) the additive cos a vanishes after the integration over angles and 
the vertex corrections vanish. 

One should note that when the transfer integral is not much greater than the cyclotron energy 
(4:Trt/huJc ~ 1) one cannot use the expansions ( |l9| ) of the Bessel's functions and should apply the 
formula ( |l8|) for the beatings of the conductivity oscillations. This slightly changes the phase shift 
and makes the beatings not periodic. This should be taken into account in the extraction of the 
transfer integral from the beating frequency. 



4 The slow oscillations 

The slow oscillations come from the interference of the oscillations of different quantities in (|Tci|). 
Hence, one should consider the entanglement of all oscillating quantities and find the non-oscillating 
term in the lowest order on damping factors. For example, the Dingle factor should no more be 
considered as constant because the entanglement of the oscillations of the Dingle factor with other 
oscillating quantities produce the term of the same order as the interference of the oscillations of 
the mean square velocity and the scattering time. The errors in the self energy now become more 
important, so, the subsequent analysis may not be quantitatively very accurate. Nevertheless, it 
describes the main features of the slow oscillations of the conductivity. 
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Up to the square of the Dingle factor in Eq. (uM we get 



e 2 N LL {ttdf 
8irkBTDh 2 LU c 



de (-n' F (e)) x 



f / huj c j /4ttA {2ne\ R 
y \ irt 1 \fuv c J \fuv c J 



(e)) x 



1 + 2J 



h/jJr 



f2ire 

\flU! c 



R D (e) + 2 J 



2ne 



R D (e) 



2nk B T D ( 4nt 



t 



Ji ( ) cos ( ^ ) Rd(c) 



ftw. 



(23) 



where i?_o(e) is now given by f|, 1 and [l]) 

2n 2 k B T D 



R D {e) 



exp 



llUJ r 



l-2Jn 



4?rf \ / 2vre 
hjj c J \hu> c 

D H ? 2 J) T COS 



R D 



Combining these two expressions we get: 

e 2 N LL (Atdf 



8nk B TDfi 2 i 



de(-n' F {e)) x (Si + S 2 ) 



where 



Si = 1 + 2 cos 



Jo 



4ttA _&Jc/ 27r 2 fc B r g \ 7 / 4vrf 



2?rf 



(24) 



(25) 



(26) 



describes the constant term and the main term of the usual SdH oscillations (they are the same as 
in (£§)), and 



S2 = R 2 d^Jq 



Ant 



Jo 



/Ant 
\huj c 



i + 2n 2 k B T D \ _ hw C j /4ttA / + 2ir 2 k B T D + f 2n 2 k B T D 



f}hJ r 



2nt \ fuv 



huj r 



(27) 



describes the slow oscillations. 

We shall consider the case At 3> hw c where our analysis based on the self energy in the form 
(||) has better accuracy and one can use the large argument expansions (|l9|) of Bessel functions. 
The formula (E7t) then simplifies: 



S 2 



D~2 



TT 2 2t 



2n 2 k B T D 



1 + a 2 s cos 



87rf 7T 

huZ ~ 2 



where 



as 



hu c ( (2Tr 2 k B T D /hu 



2nt \ l + 2ir 2 k B T D /hu 



and <p s = arctan (as) 



(28) 



(29) 
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We now have to consider the temperature and other smearing factors. But the expression (|2^) does 
not depend on electron energy (unlike the rapidly oscillating term (|2^) responsible for the SdH 
oscillations) and, hence, all types of smearing coming from the averaging over the electron energy 
do not affect the slow oscillations (they produce a factor 1). Hence, although the slow oscillations 
have a factor R 2 D squared, they do not have the temperature damping factor Rt entering the SdH 
oscillations. Therefore, the slow oscillations may be much stronger than the SdH oscillations. This 
was indeed observed at high enough temperatu re , [l5| , . 

Now collecting the formulas (|2(J), J2q) and ( |2q ) and performing the integration over e in ( p5|) 
one obtains up to the second order on damping factors: 



e 2 N LL W 2 d 2 I x , j2hw c u | — ( 'Ant 



2 cos — - \ K-^Vl + a 2 cos — +</)) R D R 



it 



8nkBT]j h 2 Lu c I \hw c J V tt 2 W \hui c 4 

Rl \ (30) 



frhj c ( , 2ir 2 k B T D 



77 2 2t V fiw. 



1 + a/1 



/ 8-Kt 7T 

\hoj c 2 



The formulas ( p0| ) and (|29|) make some predictions about the slow oscillations that can be 
compared with the experimental data. 

1) . The frequency F$ of the slow oscillations (of the conductivity as a function of the inverse 
magnetic field) is two times larger than the beating frequency. So, the distance between the two 
nearest nodes of the beatings should be equal to the period of the slow oscillations. 

2) . The amplitude of the slow oscillations is temperature independent. Of cause, if one consider 
the next terms in the expansion of the damping factor the temperature dependence will appear 
(for example, from the oscillations of the chemical potential). But these corrections are small. At 
T > ksTu the amplitude of the slow oscillations is larger than the amplitude of the SdH oscillations 
while at T < ksTo the slow oscillations are hardly distinguishable behind the SdH oscillations. 



3). The phase shift of the slow oscillations is given by (28) and (p9[). It is equal to 4> s /2 w <f>/2 
(the factor 1/2 is because the frequency of slow oscillations is two times greater than the beating 
frequency). Hence the peaks of the slow oscillations must be between the nearest peaks of the beats 
in the magnetization and conductivity. But this phase is quite sensitive to the exact expression 
of the self energy. Since we have used the approximate expression for ImS R (m, e) (formula ||) the 
substantial error in the estimate of <p s is possible. For example, the long range crystal imperfections 
substantially damp the oscillations of the DoS that so that the oscillations of the Dingle factor 
are much weaker than in formula (|24|). This long range electron scattering should be taken into 
account when the Im£^(m, e) is calculated. 

The temperature independent harmonic damping factor comes not only from the impurity 
scattering but also from the macroscopic inhomogeneities of a sample. For example, the long- 
range background potential makes the electron kinetic energy at the Fermi level slightly different 
over the sample that produces the additional temperature independent harmonic damping. The 
corresponding part of the Dingle temperature should not be included into (^|) and to the phase 
shifts ( |2l| ) and (p9|). This smearing like the temperature smearing does not affect the phase shift 
and the amplitude of the slow oscillations but causes an additional damping of the SdH or dHvA 
oscillations. To make a quantitative estimate of the amplitude of slow oscillations one can separate 
the impurity part of the Dingle temperature from the comparison of the phase-shift and then apply 



the formulas (|25|), ( |26| ) and ( |28[ ) without adjustable parameters to compare the amplitudes of the 
oscillations. 
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To summarize, we have calculated the main terms in the conductivity a zz in the layered com- 
pound in normal magnetic field. The obtained formula ([30|) describes the slow oscillations and the 
phase shift of the beatings. But the limitations of the calculation due to zero tilting angle and the 
approximate form of the self energy may lead to some quantitative deviations in the phase and the 
amplitude of the slow oscillations. 

We are thankful to A.M. Dyugaev and I. Vagner for the encouragement and stimulating dis- 
cussions. The work was supported by the EU ICN contract HPRI-CT-1999-40013, and grants 
DFG-RFBR No. 436 RUS 113/592 and RFBR No. 00-02-17729a. 



A Appendix 

A.l Transformation of the sums over LLs to the sums over harmonics 

To transform the sums over LL number into the harmonic sums we shall apply the Poisson sum- 
mation formula p2| 



00 00 (-00 

H /(«) = J2 e 2mkn f(n) dn , where a £ (n - 1; n ) 

n=n k= — 00 a 

This formula is valid for arbitrary function /(n). The sum in (0) becomes 



(31) 



d . „ / , / l^ 2 



^ |«g(e,n)| = ^2 t \ 4£ 2 - ( e - tuu r \ n 



h V V V 2 

n=0 



2) ) V \HloJ \hw c 



= —huj r y dn exp 2 

^ , ,k f2Trike\ f°° , . , 
= —rku c y (— 1) exp — / ax exp {2m kx 



hio, 



X 2 



4^ dt , , k / 2nike\ ( A-Kkt\ 
£ ilHfexp — J: - 



k= — oo 

In an analogous way one can also transform the sum 



(32) 



N 



oo 1 

^Po(^) = E / ( 33 ) 
„=o \JU 2 - (e - 7kj c (n+ I)) 
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1 ~ dn exp (2mk (n - |)) 



( 2t 



1 \ -> , ,fc / 2i:ike\ f°° dx exp (2i:ikx) 

L »p ( — ) / 7= , 2 



exp 



h0J c J— 1 \ hu! c J \ tlUJ c J 
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